Disorder to chaos transition in the conductance distribution of corrugated waveguides 
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We perform a detailed numerical study of the distribution of conductances P(T) for quasi-one- 
dimensional corrugated waveguides as a function of the corrugation complexity (from rough to 
smooth). We verify the universality of P(T) in both, the diffusive ((T) > 1) and the localized ((T) <§; 
1) transport regimes. However, at the crossover regime ((T) ~ 1), we observe that P(T) evolves 
from the surface-disorder to the bulk-disorder theoretical predictions for decreasing complexity in 
the waveguide boundaries. We explain this behavior as a transition from disorder to deterministic 
chaos; since, in the limit of smooth boundaries the corrugated waveguides are, effectively, linear 
chains of chaotic cavities. 

PACS numbers: 73.23.-b, 71.30.+h, 42.25.Dd 



I. INTRODUCTION 

In studies of wave propagation through disordered 
wires two setups are mostly considered: bulk-disordered 
and surface-disordered waveguides. In both cases it is 
possible to discern between diffusive (metallic) and lo- 
calized (insulating) transport regimes by comparing the 
wire length L with the mean-free path £ and the local- 
ization length £; i.e., wave diffusion takes place when 
£ <C L <C £ while localization is observed for £ <C L. 

From the analytical point of view both disorder setups 
have been successfully approached. On the one hand, 
transport through bulk-disordered waveguides is well de- 
scribed by the Fokker-Planck approach of Dorokhov, 
Mello, Percyra, and Kumar (DMPK) U@|; and by the 
field-theoretic approach of Efctov and Larkin, which 
leads to a supersymmetric nonlinear a model In 
fact, these two approaches were shown to be equivalent in 
Ref. pi| • In addition, in Ref. [ll| , the distribution of con- 
ductances P(T) in the full diffusive-to-localized crossover 
was derived, in the frame of the supersymmetric ap- 
proach, for waveguides with broken time-reversal invari- 
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On the other hand, transport through surface- 
disordered wires has been properly characterized by 
the Fokker-Planck approach developed by Froufe-Perez, 
Yepez, Mello, and Saenz (FYMS) [13j . Other analytical 
approaches to transport throu gh s urface-disordered wires 
are also available in Refs. flil Il5l|. 

Furthermore, for both setups, P(T) evolves from a 
Gaussian shape (deep in the diffusive regime) to a log- 
normal shape (deep in the localized regime). However, 
at the crossover between diffusive and localized regimes, 
the form of P(T) is highly non-trivial @, [IH, Il6l - [23j and 
importantly depends on the type of disorder (bulk or 
surface). Moreover, the DMPK approach and the FYMS 
approach provide accurate predictions for P(T) at the 
crossover regime for the corresponding setups of disorder 

urn EM. 

In this paper we numerically study P(T) for a model 
of quasi-one-dimensional surface-disordered waveguides 
with tunable corrugation complexity: from rough to 
smooth. Here, we concentrate on waveguides with time- 



reversal invariance. We define the corrugated surface of 
our disordered wire as a sum of harmonics with random 
amplitudes. In the rough limit (large number of harmon- 
ics) the waveguide effectively shows surface disorder; so 
that, it is equivalent to the step-like corrugated waveg- 
uide model used in Refs. [16Hl9l l. On the other hand, 
in the smooth limit (few harmonics) the waveguide can 
be considered as a linear chain of attached chaotic cav- 
ities. Interesting enough, the transport properties of a 
waveguide constructed as a linear chain of chaotic cavi- 
ties l9l.[23 . [25| are equivalent to those of a bulk-disordered 
wire [IOJ. Then, by decreasing the corrugation complex- 
ity of our surface corrugated waveguide we expect to ob- 
serve, at the diffusive-to-localized transition regime, a 
transition in the form of P(T) from the surface-disorder 
FYMS to the bulk-disorder DMPK predictions. 

The organization of this paper is as follows. In the 
next section we define the waveguide model we use as 
well as the scattering setup. In Sec. Ill, by extracting £ 
and £ from curves of the average resistance and average 
logarithm of conductance as a function of L, respectively, 
we define the diffusive and localized transport regimes 
for our corrugated waveguides. Then we study in detail 
the distribution of conductances P(T) as a function of 
the corrugation complexity in both regimes, diffusive and 
localized, as well as at the crossover regime. Finally, Sec. 
IV is left for conclusions. 



II. MODEL 



The model we shall use in our study is a waveguide 
formed by attaching L two-dimensional cavities. Each 
cavity of length L x is defined by two hard walls: one 
flat at y = and a corrugated one given by y = d + 
e/(x). Here d is the average width of the cavity and 
e is the corrugation amplitude. Since we are interested 
in studying the transport properties of waveguides as a 
function of the complexity of the corrugated boundary 
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FIG. 1: Examples of cavities used to form corrugated waveg- 
uides. Here we show one realization of cavities with Nt = 20 
and Nt = 3. In our study we go from Nt = 18 to Nt = 2. 



The experimental realization of a scattering setup sim- 
ilar to ours has been recently reported in Ref . (32[ . How- 
ever, we notice that in numerical simulations of trans- 
port through surface-disordered wires, step- like corru- 
gated waveguides are more often used (l6Hl9l . l33U35j [. 
among others [HI, |36| . We also note that here we con- 
centrate on the case of small number of open modes, 
M = [2,9]; the case of M ^> 1 has been recently ad- 
dressed in Ref. 35]. 



we choose 

f( x ) = f^A n cos(?p x \ , (1) 

n— 1 \ x / 

where A n are random numbers drawn from a flat distri- 
bution in the range [—1,1]. This form for f(x) allows 
us to choose the desired degree of complexity of the cor- 
rugated waveguide boundary: from rough, Nt ~ 20, to 
smooth, Nt ~ 1. Sec Fig. [T] It is important to stress 
that once the parameters L x , d, e, and Nt are fixed we 
randomly generate (through different values of A n ) L 
different cavities that we attach to form a non-periodic 
waveguide 0, [27} . To this end the minimal Nt we use 
is 2. 

We remark that depending on the values of the param- 
eters (L x , d, e, and Nt), the classical (or ray) dynam- 
ics in each of the cavities can yield mixed or full chaos. 
However, here we consider only the case of full chaos 
[13: edN T /L 2 x > 0.01. Then, below we use L x = 2tt, 
d = L x /2, and e = L x /20. This set of parameters pro- 
duces full chaos for any Nt > 1. Finally, note that all 
lengths, here and below (including £, £ and L), are given 
in units of L x . 

We open the waveguide of length L defined above by 
attaching two semi-infinite collincar fiat leads of width d 
to its left and right ends. The leads support plane waves 
with energy E = (?i 2 /2m)[fc^ + (mir/d) 2 ], where k m and 
mir/d are, respectively, the longitudinal and transver- 
sal components of the total wave vector K = v / 2mB/S. 
Then, using finite element methods (see e.g. [Ij|[3(J) we 
compute the scattering matrix, S'-matrix, which has the 
form 

»-(;?)• 

where t, t 1 , r, and r' are M x M transmission and re- 
flection matrices. Here, M is the highest mode given by 
the largest m beyond which the longitudinal wave vector 
k m = [2mE/h 2 — (mir/d) 2 ] 1 / 2 becomes complex. Then, 
once the S'-matrix is known we calculate the dimension- 
less conductance from [3l[ 

T = Tr(tt f ) . (2) 

With this definition, the conductance can take values in 
the interval [0, M]. 



III. RESULTS 
A. Diffusive and localized regimes 

In order to identify the diffusive and localized trans- 
port regimes in our corrugated waveguides, in Fig. [2] we 
plot the average resistance (1/T) and the average loga- 
rithm of conductance (InT) as a function of L for waveg- 
uides with N T = 18, 8, 3, and 2. 

For disordered wires, it is well established that (i) for 
relatively short wire lengths the resistance increases lin- 
early with L as [HI, I36l - l38j 

1 \ _ 1 L 

T/M + M£ ] () 

while (ii) for relatively large wire lengths the conductance 
decays exponentially with L in the form [ltl HH, H^, [4f| 

(lnT)oc-|. (4) 

We extract I and £ from Fig. [2] by performing fittings 
of the data with Eqs. © and ([J}, respectively, see red 
dashed lines. In Fig. [5] we also indicate, with blue dot- 
dashed vertical lines, the positions of the obtained £ and 
£. From this figure, it is clear that both £ and £ decrease 
by increasing Nt] cquivalcntly, the mean- free path and 
the localization length increase when the disorder de- 
creases. Therefore, in Fig. [2] we label as (II) and (III) 
the diffusive and localized regimes, respectively. Addi- 
tionally, we identify with (I) the quasi-ballistic regime, 
L < £, which we will not explore here. 

It is also interesting to mention that we found that 
relation £ rj Ml (see for example [l?], Ell) works well 
only when the surface corrugation is complex enough; i.e. 
for Nt > 7. See upper panels of Fig. [2] When Nt — > 1 
we observe that £ < Ml. See lower panels of Fig. [5J In 
any case the diffusive regime is clearly discernible in our 
calculations for all the values of M and Nt we used here. 

Note that to construct Fig. [2] we have used waveguides 
supporting M — 3 open modes (here and below, when 
showing results for different values of M we always fix 
the energy such that K = n(M + l/2)/d). We obtained 
similar plots for other values of M. However, by increas- 
ing M the values of £ and £ decrease and, as a conse- 
quence, the regimes (I) and (II) become narrower. For 
M = 9, the highest M we explored, the quasi-ballistic 
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FIG. 2: (Color online) Average resistance (1/T) (upper panels) and average logarithm of conductance (lnT) (lower panels) as 
a function of the length L of corrugated waveguides with (a) Nt = 18, (b) 8, (c) 3, and (d) 2 (black full lines). To average, 
10 4 waveguide realizations were used. The waveguides support M = 3 open modes. Red dashed lines are best fittings of 
(1/T) [(lnT)] with Eq. © [Eq. Q}] for small [large] L. The extracted (£,£) are approximately equal to (a) (44.77,134.38), 
(b) (55.88, 179.28), (c) (94.97, 265.78), and (d) (136.88, 346.61). Blue dot-dashed vertical lines indicate the positions of i and £ 
which delimit the transport regimes: quasi-ballistic (I), diffusive (II), and localized (III). Red vertical lines mark the value of 
ML 



regime is hardly visible in the scale of Fig. [2] In the fol- 
lowing figures we will use values of M different from 3 to 
emphasize that our results do not depend on the num- 
ber of open modes, once the transport regimes are well 
determined. 

Then, in Fig. [3l we verify the predictions for the con- 
ductance probability distribution function in the diffu- 
sive and localized regimes. As clearly shown in this fig- 
ure, where M = 9 open modes were considered, P(T) 
has a well defined Gaussian shape in the diffusive regime 
[Fig. Elja)], while P(lnT) has the log-normal form in the 
localized regime [Fig. [3fb)]. We know that the results 
reported in Fig. [3] are already expected and thus it may 
seem unnecessary to shown them. However, we decided 
to present them in order to stress that the transport prop- 
erties, exemplified here by the form of P(T), of rough and 
smooth corrugated waveguides are in fact similar except 
at the crossover regime, as will be shown below. 



B. Crossover regime 

As already anticipated, we have found that the most 
interesting result for the distribution of conductances ap- 



pears at the crossover regime, where (T) ~ 1. Here, 
we observe that the shape of P(T) does depend on the 
waveguide corrugation complexity. To show this, in Fig. [4] 
we plot P(T) for corrugated waveguides with Nt — 18, 8, 
3, and 2 (columns) for some values of (T) (for comparison 
purpos es in Fig. [4] we chose the same values of (T) used in 
[la . fla - flSll ] ) . As a reference, we also include the surface- 
disorder FYMS and the bulk-disorder DMPK predictions 
for P(T). In each panel we plot three histograms corre- 
sponding to waveguides supporting M = 5,7, and 8 open 
modes. Note that for large Nt, in our case Nt = 18 |41| . 
the shapes of the numerically obtained P(T) are well de- 
scribed by the FYMS prediction for surface-disordered 
waveguides, as expected (see panels in the left-most col- 
umn of Fig. |4]). However, once the waveguide corruga- 
tion complexity is decreased, important deviations ap- 
pear. Moreover, when Nt = 2, P(T) fully coincides 
with the DMPK prediction for bulk-disordered waveg- 
uides (see panels in the right- most column of Fig. |4j. 
This fact is more evident for (T) = 1 and 4/5 where the 
differences between FYMS and DMPK predictions are 
easily distinguishable. 

So, we observe an effective and smooth evolution in 
the form of P(T), from the surface-disorder FYMS to 
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FIG. 3: (Color online) Conductance probability distributions 

(a) P(T) and (b) P(lnT) for corrugated waveguides with 
Nt = 18, 8, 3, and 2 (from left to right) in the (a) diffusive 
and (b) localized regimes. 10 4 values of T where used to con- 
struct each histogram. The waveguides support M = 9 open 
modes. Dashed lines are (a) Gaussian and (b) log-normal dis- 
tribution functions characterized by the values of (a) (T) and 

(b) (In T) reported in the corresponding panels. 



the bulk-disorder DMPK predictions, as a function of 
the (decreasing) waveguide corrugation complexity Nt- 
We understand this result in the following simple way. 
In the limit Nt 1 our waveguide can be considered 
as a linear chain of coupled chaotic cavities, each one 
defined by the cosine billiard Then, the model of 

quantum dots in series of Iida, Weidenmuller, and Zuk 
applies [HI ■ Moreover, this model reduces to a super- 
symmetric nonlinear a model 0-Q which turns out to be 
equivalent to the DMPK approach That is, while 

FYMS describes well the limit Nt ^ 1 of our corrugated 
waveguide, DMPK should describe the limit Nt — > 1; as 
we in fact observe in Fig. 2) Therefore, our waveguide 
model shows a disorder-to-chaos transition in the shape 



of the conductance distribution. 

Fin ally , we want to add that the parameters used in 
Refs. [l6Hl9l |. translated to our symbols, are [43[: d/e = 
7 - 13.25 (here we used d/e = 10) and N T « 7. Then, 
since we are observing well developed surface-disorder 
FYMS transport properties in our corrugated waveguides 
for Nt > 8, there is no contradiction between our results 
and those presented in Refs. (l6l - fl9| . Moreover, we think 
that it would be interesting to explore the limit where the 
step length in the step-like corrugated waveguide model 
of Refs. |16l - fl9l |. becomes of the order of the waveguide 
width; which is somehow equivalent to the limit Nt — > 1 
in our waveguide model. 



IV. CONCLUSIONS 

We have studied the distribution of conductances P(T) 
for a quasi-onc-dimcnsional corrugated waveguides with 
tunable corrugation complexity: from rough, Nt = 18, 
to smooth, Nt = 2. We verified that both, the mean free 
path and the localization length decrease for increasing 
N T - Also, we confirmed that P(T) and P(lnT) have the 
Gaussian and log-normal shapes in the diffusive ((T) > 1) 
and localized ((T) -C 1) transport regimes, respectively. 

At the crossover between the diffusive and the localized 
regime, (T) ~ 1, we reported that P(T) monotonously 
evolves from the surface-disorder FYMS (when Nt = 18) 
to the bulk-disorder DMPK (when Nt = 2) predictions 
for decreasing Nt- We understood this behavior as a 
consequence of the underlying deterministic dynamical 
chaos; since the waveguides having smooth boundaries 
(i.e., when Nt — >• 1) are effectively linear chains of at- 
tached chaotic cavities. 

We believe that our results, as well as our model of cor- 
rugated wires with tunable corrugation complexity, may 
stimulate further analytical and numerical studies on the 
transport properties at the diffusive-to-localized transi- 
tion regime. 
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